
%% initialize
clear all;
close all;

rng(328120)

Nsim = 500;

loadstarting = 0;
loadestimates = 0;

%% read in data

% choices

load '../../01_data/05_temp/mixed_logit_rc6_totVA_DISAD_Ndis'


rawdata = csvread('../../01_data/02_basedata/estimationdata/All_T_J/All_T_J_choices_buff_0_0_lf_tbar_Ndis_C.csv',1,0);


IJVAma          = rawdata(:,22);
IJVAmean        = rawdata(:,23);

IJVAm1          = rawdata(:,end-3);
IJVAm2          = rawdata(:,end-2);


clearvars rawdata;

rawdata = csvread('../../01_data/02_basedata/estimationdata/All_T_J/All_T_J_choices_buff_0_0_lf_tbar_Ntype_C.csv',1,0);
N1              = rawdata(:,3);
N2              = rawdata(:,4);
IJfracDisadv = N2./(N1+N2);
IJclassSize = N1+N2;

clearvars rawdata;



 save -v7.3 '../../01_data/05_temp/mixed_logit_rc6_totVA_DISAD_Ndis_Calt'
 
